# simulate a dataset with continuous data
y = rdag(1000, 15, 0.3)
tru = y$G
x = y$x
mod = pc.con(x)
mod$runtime
b = pc.or(mod)
plota(tru)
dev.new()
plota(b$G)
dev.off()
plota( dag2eg(tru) ) ## essential graph
dev.new()
plota(b$G)
Run the code above in your browser using DataLab